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We present results for the static potential in the confinement phase of the SU(2) Higgs model on the lattice. 
Introducing a suitable matrix correlation function we observe string breaking at a distance « 1.8ro, where the 
length scale ro has a value ro ~ 0.5 fm in QCD. The method presented here may lead the way towards a treatment 
of string breaking in QCD. 



1. Introduction 

In pure non-Abelian gauge theories there is 
a linear confinement potential between a static 
source anti-source pair in the fundamental repre- 
sentation of the gauge group at large distances. 
The static potential can be efficiently computed 
by means of Wilson loops. C. Michael JIJ] stud- 
ied the potential between static adjoint sources 
in the pure SU(2) gauge theory: the gauge fields 
screen the sources forming color-neutral objects 
called gluelumps. At large distances the two 
gluelump state dominates the ground state of the 
system. A similar situation is expected when 
these gauge theories are coupled to matter fields 
in the fundamental representation. The dynam- 
ical matter fields form a bound state with the 
static fundamental source, a color-neutral static 
"meson". One does expect that the potential at 
large distances is better interpreted as the poten- 
tial between two static mesons and flattens turn- 
ing asymptotically into a Yukawa form. So far, 
this expectation could not be verified by Monte 
Carlo simulations. In particular, in recent at- 
tempts in QCD with two flavors of dynamical 
quarks this string breaking effect was not visible 
(2j||] . The distance r b around which the potential 
should start flattening off, could be estimated in 
the quenched approximation pi: 



r b «2.7r (in QCD) 



(1) 



where ro is the reference scale defined in Q . This 
is in agreement with an estimate from full QCD 
simulations 0]. The same behavior of the po- 



tential is, of course, expected in the non-Abelian 
Higgs model in the confinement phase. While the 
potential in the large distance range could not be 
calculated in early simulations with gauge group 
SU(2), they yielded some qualitative evidence for 
screening of the potential 1^,0. In our work, we 
compute the potential in the Higgs model and ob- 
serve string breaking. A recent study in the three- 
dimensional SU(2) Higgs model || has reached 
very similar conclusions to what we find in four 
dimensions. 

2. Calculation of the potential 

For our investigation of the static potential in 
the SU(2) Higgs model we choose in the conven- 
tional bare parameter space Q the point f3 — 2.2, 
k = 0.274 and A = 0.5. It lies in the confinement 
phase of the model, fairly close to the phase tran- 
sition, where the properties are similar to QCD. 
The lattice resolution is of roughly the same size 
as the one used in the QCD-studies: we obtain 
r /a = 2.78 ±0.04. 

The general strategy for determining the static 
potential has first been applied in 0. The cor- 
relation functions that we used are schematically 
illustrated in Fig. [j]. We measure a symmetric 
matrix correlation function Cij (r, t) , where the 
indices i and j refer to the space-like parts of the 
correlation functions. These consist of a Wilson 
line or two Higgs fields. For fixed spacial separa- 
tion r, the potential V(r) is extracted solving the 
generalized eigenvalue problem pQ] : 



C(t)v a (t,t ) 



X a (t,to)C(t )v a (t,to), (2) 
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with X a > Aq,+i. The ground state energy V(r) = 
Vo(r) and the excited states energies Vi(r), ... are 
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Figure 1. The correlation functions used to deter- 
mine the static potential. The lines represent the 
Wilson lines, the filled circles the Higgs field. 



then given by 

aV a (r) = ln(A Q (< - a,t Q )/\ a (t, t Q )) 



-0 



-(V N (r)-V a (r))t 



(3) 



Here, N is the rank of the matrix C. In order 
to suppress the correction term in eq. (^) we use 
smeared gauge and Higgs fields at different smear- 
ing levels. In this way the rank of C is increased 
to N = 7. For more details see J]~T| ] . 

3. String breaking 

Our numerical results were obtained on a 20 
lattice with periodic boundary conditions. We 
computed all correlation functions up to r max = 
imax = 8a ~ 3 ro on 4240 field configurations. 
Statistical errors were reduced by replacing - 
wherever possible - the time-like links by the 1- 
link integral. 

We extracted the mass /i of a static meson using 
the variational method of |l(J from a correlation 
function with one straight time-like Wilson line 
and smeared Higgs fields at the ends. In Fig. 2 
we compare two different smearing procedures for 
the Higgs field. The smearing operator iSi (trian- 
gles) is defined as: 

5i$(i) - 9(x) + ]T U{x,y)*{y), (4) 

|*-y|=« 

xo=yo 

where is the complex Higgs field and U(x, y) 
is the gauge field link connecting y with x. The 
smearing operator S2 (circles) is defined as: 

S 2 ${x) = V{V^>(x)+V U{x,y)<$>{y) 

|x- B | = V2a 
^0 = V0 
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Figure 2. Comparison of the extraction of the mass /i 
of a static meson using the smearing operators defined 
in eqs. (§) and (§) . 



-V ]T U(x,z)$(z)}, 



(5) 



\=V3a 



where — $/v $t$ and U(x,y) represents 
the average over the shortest link connections be- 
tween y and x. The application of the smear- 
ing operator is iterated obtaining a sequence 
c[)(ro) ^ _ s m <&(x) with which the correlation 
function is evaluated. Using £2 the contribu- 
tions from the excited states are much more sup- 
pressed: at t = 7a we read off afi = 0.7001 ± 
0.0014 which agrees fully with t = 6a. 

We computed the potential in units of ro . Con- 
sidering in particular the combination V(r) — 2[i, 
one has a quantity free of divergent self energy 
contributions. It is shown in Fig. y. The expected 
string breaking is clearly observed for distances 
r > rb ~ 1.8 ro- Around r ~ r^, the potential 
changes rapidly from an almost linear rise to an 
almost constant behavior. 

We want to point out that if one considers only 
the sub-block of the matrix correlation function 
corresponding to the (smeared) Wilson loops, the 
potential estimates have large correction terms at 
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Figure 3. The static potential in units of ro. The 
asymptotic value 2p has been subtracted to obtain 
a quantity free of self energy contributions which di- 
verge like -. 



long distances. One might then extract a poten- 
tial which is too high. This might explain why 
string breaking has not been seen in QCD, yet. 
This observation is confirmed by the overlap of 
the variationally determined ground state, char- 
acterized by Vq in eq. (^), with the true ground 
state of the Hamiltonian. The overlap can be 
computed from the projected correlation function 



Q(t) = v^C(t)v 



oj n e 



-V n (r)t 



(0) 



with normalization 51(0) = 1. Here, n labels the 
states in the sector of the Hilbert space with 2 
static sources. "The overlap " is an abbreviation 
commonly used to denote ujq. Considering the 
full matrix C the overlap exceeds about 50% for 
all distances, restricting the matrix to the Wilson 
sub-block we find upper limits for the overlaps at 
r>r h of 5% 0. 

4. Conclusions 

We have introduced a method to compute the 
static potential at all relevant distances in gauge 



theories with scalar matter fields. We demon- 
strated that it can be applied successfully in the 
SU(2) Higgs model with parameters chosen to re- 
semble the situation in QCD. It is then interest- 
ing to follow a line of constant physics towards 
smaller lattice spacings in order to check for cut- 
off effects and to be able to resolve the interesting 
transition region in the potential. 

From the matrix correlation function one can 
also determine excited state energies. A precise 
determination of the excited potential at all dis- 
tances needs more statistics. One expects that 
the transition region of the potential can be de- 
scribed phenomenologically by a level crossing (as 
function of r) of the "two meson state" and the 
"string state" We are planning to investi- 

gate this in more detail. So far, we can only say 
that for r w rb the two levels V\ (r) and V(r) are 
close. 

Of course, it is of considerable interest to apply 
this method to QCD with dynamical fermions. 
The only possible difficulty is expected to be one 
of statistical accuracy. The proper correlation 
functions can be constructed along the lines of 
ref. 0,1,0. 

REFERENCES 

1. C. Michael, Nucl. Phys. B (Proc. Suppl.) 26 
(1992) 417. 

2. R. Burkhalter for CP-PACS Collaboration, 
these proceedings. 

3. C.R. Allton et al., (1998), |hep-lat/9808016| . 

4. R. Sommer, Nu cl. Phys. B411 (1994) 839, 
hep-lat/9310022 . 

5. R. Sommer Phys. Rept. 275 (1996) 1, [hep 



lat/9401037 



6. H.G. Evertz et al., Phys. Lett. 175B (1986) 
335. 

7. W. Bock et al., Z. Phys. C45 (1990) 597. 
O. Philipsen and H. Wittig, (1998), [hep" 



lat/980702C 



9. I. Montvay, Nucl. Phys. B269 (1986) 170. 

10. M. Liischer and U. Wolff, Nucl. Phys. B339 
(1990) 222. 

11. F. Kncchtli and R. Sommer, (1998), [hep 
lat/9807022] 



12. IT. Drummond, (1998), |hep-lat/9805012| . 



